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ABSTRACT 

Hot (explosive) hydrogen burning or the Rapid Proton Capture Process (rp-process) 
occurs in a number of astrophysical environments. Novae and X-ray bursts are the most 
prominent ones, but accretion disks around black holes and other sites are candidates 
as well. The expensive and often multidimensional hydro calculations for such events 
require an accurate prediction of the thermonuclear energy generation, while avoiding 
full nucleosynthesis network calculations. 

In the present investigation we present an approximation scheme which leads to an 
accuracy of more than 15 per cent for the energy generation in hot hydrogen burning 
from 10 8 -1.5 x 10 9 K, which covers the whole range of all presently known astrophysical 
sites. It is based on the concept of slowly varying hydrogen and helium abundances 
and assumes a kind of local steady flow by requiring that all reactions entering and 
leaving a nucleus add up to a zero flux. This scheme can adapt itself automatically and 
covers situations at low temperatures, characterized by a steady flow of reactions, as 
well as high temperature regimes where a (p, 7 )~( 7 , p)-equilibrium is established, while 
/3 + -decays or (a,p)-reactions feed the population of the next isotonic line of nuclei. 

In addition to a gain of a factor of 15 in computational speed over a full network 
calculation, and an energy generation accurate to more than 15 per cent, this scheme 
also allows to predict correctly individual isotopic abundances. Thus, it delivers all 
features of a full network at a highly reduced cost and can easily be implemented in 
hydro calculations. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — stars: novae - 
X-rays: bursts 
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1. INTRODUCTION 

Close binary stellar systems can exchange 
mass, when at least one of the stars fills its 
Roche Lobe. After a critical mass AM of 
unburned transferred matter is accumulated 
on the surface of the accreting star, igni¬ 
tion sets in, typically under degenerate con¬ 
ditions when the accreting object is a white 
dwarf or neutron star. Degenerate condi¬ 
tions for which the pressure is not a function 
of temperature, prevent temperature adjust- 
mq nt via pressure increase and expansion and 


cai fsp a thermonuclear r.im-nTm—rrni I exploswe- 

bu: 


ning. — Foi while dwarfs, a fever of KU¬ 

RD 4 M 0 forms, before pycnonuclear ignition 
of hydrogen burning sets in (e.g. [Sugimoto 
fe Nomoto 198Q|; (Stanfield ct al. 1993 ; Coc 


t ftl, 1995| ) and causes a nova event , (explo¬ 


siv e Il-burning on white dwarfs) with maxi- 

mum temperatures of 


;(2-3)xl0 8 

^46 1 n47 


K and a 


total energy release of 10 4 -10 4 ' ergs. The 
burning takes 100-1000 s before the partially 
burned hydrogen envelope is ejected. 

X-ray bursts (for an observational overview 
see Lewim van Paradijs, fe Taam| 1993) in¬ 
volve accreting neutron stars with an un¬ 
burned, hydrogen-rich surface layer. The 
critical size of the hydrogen layer before ig¬ 


nition is as small as 10 


-12 


M, 


©■ 


Temper¬ 


atures of the order (1-2) xlO 9 K and den¬ 
sities p«H0 6 -10' g cm -3 are attained (see 
e.g. |Wallacc fc Wooslcy 1981[ |Aya,sli fe Joss 
1982; [Hanawa, Sugimoto, fe Hashimotoj 1983; 
Woosley, Taarn, fe Weaver| 1986; [Taarn et~ah 


1993; [Th am. Wooslcy. k, Lamb| 1996). This 
explosive burning with rise times of about 1- 
10 s leads to the release of 10 39 -10 4 ° ergs. 
Many of the observed features and general 
characteristics are understood, however, there 
is still a lack of a quantitative understand¬ 
ing of the detailed observational data. An¬ 


other aspect is that the explosion energies are 
smaller than the gravitational binding energy 
of the accreted hydrogen envelope. An evenly 
distributed explosion energy would not un¬ 
bind and eject matter. It remains to be seen 
whether an interesting amount of matter can 
escape the neutron star. Super-Eddington X- 
ray bursts ( [Taam ct al. 1996|) are the best 
candidates for this behavior. 

A description of hot (explosive) hydrogen 
burning has been given by |Arnould ct al.| 
(1980), [Wallace fe Wooslcyl (1981), |Ayasli & 
Joss| (1982), [Hanawa et al.| (1983), |Wicschcr| 
ct al.| (1986), |van Wormcr ct al] (1994), 
[rhielemann et aT| (1994) [see also [Bichlc| 
(1991), |BichIc| (1994), |Cannon et al.| (1992), 
and |(Jannon| (1993) for Thorne-Zytkow ob¬ 
jects, but see the recent results by [Fryer ,| 
Benz, fe Herant| (1996) which puts doubts 


on their existence]. The burning is described 
by proton captures, /D-decays, and possibly 
alpha-induced reactions on unstable proton- 
rich nuclei, usually referred to as the rp- 
process (Rapid Proton Capture Process). Cross 
sections can either be obtained from the best 
available application of present experimen¬ 
tal knowledge, e.g. the determination of 
resonance properties from mirror nuclei and 
transfer and/or charge exchange reaction stud¬ 
ies, or actual cross section measurements, like 
for 13 N(p, 7) 14 0, the first reaction cross sec¬ 
tion analysed with a radioactive ion beam fa¬ 
cility (see IChampagnc fc Wicschcr| 1992 and 
references therein). There exist two major 
motivations in nuclear astrophysics (a), to 
understand the reaction flow to a necessary 
degree, in order to predict the correct en¬ 
ergy generation required for hydrodynamic, 
astrophysical studies, and (b), to predict a 
detailed isotopic composition which helps to 
understand the contribution of the process in 
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question to nucleosynthesis in general. Our 
main motivation in this paper is (a), to pro¬ 
vide a fast energy generation network as a 
tool for nova and X-ray burst studies. This 
might be underlined by the fact that no¬ 
vae and X-ray bursts seem not to be ma¬ 
jor contributors to nucleosynthesis due to the 
small ejected masses involved or the question 
whether gravitational binding can be over¬ 
come at all. [However, novae can be impor¬ 
tant for nuclei like 'Li, 15 N, 22 Na, 26 A1, and 
even Si and S; and super-Eddington X-ray 
bursts ([Taam et al.| 1996) will be able to eject 
some matter, probably containing some light 
p-process elements.] It will turn out at the 
end that our efficient approximation and en¬ 
ergy generation scheme can also be used to 
predict abundances accurately. The nucle¬ 
osynthesis for nuclei above Kr will, however, 
be discussed in a second paper (|Schatz| et al. 
1996). 

In order to understand the energy gener¬ 
ation correctly, we have to be able to under¬ 
stand the main reaction fluxes. In the past we 
performed a series of rp-process studies ( |Wi- 


esc jher et al.| 1986 to Ar, and [van Worrner 
et ah 1994, extending up to Kr). The reac¬ 
tion rate predictions were based on resonance 
and direct capture contributions for proton- 
rich nuclei, making use of the most recent ex¬ 
perimental data. Several (p, 7 )-reaction rates 
below mass A=44 have been recalculated by 
Herndl et al.| (1995) in the framework of a 
shell model description for the level structure 
of the compound nucleus. A preliminary anal¬ 
ysis of the major aspects was given by |Thiclc- 


me nn 


et al.| (1994). In the following section 


H we will present this in more detail and dis¬ 
cuss the constraints which an approximation 
scheme has to fulfill in the whole range of 
temperatures occurring in explosive hydrogen 


burning environments, i.e. 10 8 -1.5 x 10 9 K. 
The approximation scheme will be presented 
in section |3|, its application and comparison 
with full network calculations in section [|. fol¬ 
lowed by concluding remarks in section |j. 

2. REACTION FLOWS IN EXPLO¬ 
SIVE HYDROGEN BURNING 

At low temperatures, the rp-process is dom¬ 
inated by cycles of two successive proton cap¬ 
tures, starting out from an even-even nucleus, 
a /? + -decay, a further proton capture (into an 
even-Z nucleus), another /3 + -decay, and a fi¬ 
nal (p, a)-reaction close to stability (similar 
to the hot CNO). The final closing of the cy¬ 
cle via a (p, a)-reaction occurs because pro¬ 
ton capture Q-values increase with increas¬ 
ing neutron number N in an isotopic chain 
for a given Z, while a-capture Q-values or - 
thresholds show a very weak dependence on 
N. Thus, a cycle closure occurs at the most 
proton-rich compound nucleus, which has a 
lower alpha than proton-threshold. These 
are the even-Z and especially well bound 
“alpha”-nuclei with isospin T Z =(N — Z)/2=0, 
like 16 0, 20 Ne, 24 Mg, 28 Si, 32 S, 36 Ar, 40 Ca. 
Therefore, (p, a)-reactions can operate on odd- 
Z targets with T*=+l/2: 15 N, 19 F, 23 Na, 
27 Al, 31 P, 35 C1, 39 K. The hydrogen burning 
cycles are connected at these nuclei, due to a 
possible competition between the (p, 7 )- and 
(p, a)-reaction. An exception to the rule is 
the Ne isotopic chain, where the transition 
happens already for the compound nucleus 
19 Ne. Thus, 18 F(p, a) 15 0 is a possible reac¬ 
tion, bypassing the OF(Ne)-cycle. Therefore, 
only the extended CN(O)-, the NeNa-, MgAl- 
, SiP-, SCl-cycles etc. exist. This is displayed 
in Figure [|. 

The progress of burning towards heavier 
nuclei depends on the leakage ratio (p, 7 )/(p, a) 
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Fig. 1.— The hot hydrogen burning cycles, Fig. 2.— Q-values of (p, 7 )-reactions between 

established at T g =0.3 and p= 10 4 g cm -3 . 20 Ne and 44 Ti. 


into the next cycle, which makes a good ex¬ 
perimental determination of these reactions 
important and is possible as we deal here with 
stable targets. On the other hand, the excita¬ 
tion energies in the corresponding compound 
nuclei are of the order 8.5-12 MeV and make 
statistical model approaches also reliable (see 
e.g. the experiments by |Iliadis ct al.| 1993, 
1994 and [Ross et ah] 1995 for 31 P(p, 7 ), and 
35 Cl(p, 7 ) and the later discussion). 

Increasing temperatures allow to overcome 
Coulomb barriers and extend the cycles to 
more proton-rich nuclei, which permit addi¬ 
tional leakage via proton captures, compet¬ 
ing with long /3-decays. The cycles open first 
at the last /3-decay before the (p, a)-reaction 
via a competing proton capture, i.e. at the 
even-Z T z =-l/2 nuclei like 23 Mg, 27 Si, 31 S, 
35 Ar, and 39 Ca, excluding 15 0 (because 16 F 
is particle unstable) and 19 Ne, which is by¬ 
passed by the stronger 18 F(p, a)-rate. These 
nuclei are only one unit away from stability 
and have small decay Q-values and long half- 
lives. The situation is illustrated in Figure 


j 2 ], which shows the possible break-out points 
and the Q-values of the proton capture reac¬ 
tions. Since the reaction proceeds towards the 
proton drip-line and into an odd-odd nucleus, 
the capture Q-values of these (p, 7 )-reactions 
are small (generally less than 2 MeV). This 
implies a too small density of excited states 
for employing statistical model cross sections 
(see discussion in the following paragraphs). 
At about 3 x10 s K essentially all such cycles 
are open and a complete rp-pattern of proton 
captures and /3-decays is established, with the 
exception of the extended CNO-cycle. The 
process time is determined by the sum of 
/3-decay and proton capture time scales, 77 
and r P)7 , along the rp-path, which does gener¬ 
ally not extend more proton-rich than T z =-1 
for even-Z nuclei. The reaction time scales 
are dominated by the T z =-l/2 even-Z nuclei 
close to stability, listed above. 

At 4x10 s K the CNO-cycle breaks open 
via 15 0(a, 7 ) 19 Ne. 14 0(ct,p) 17 F can also lead 
to a successful break-out when followed by a 
proton capture and (a, p )-reaction on 18 Ne. 
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Due to higher energies, Coulomb barriers can 
be more easily overcome by heavier projec¬ 
tiles (like alpha captures on CNO-targets), 
and proton captures can occur further away 
from stability, ffowever, in all cases the decay 
of T,= 1 even-Z nuclei (and at high densi¬ 
ties the competition with (a,p)-reactions) is 
a prominent part of the reaction path. Alpha- 
reactions are not important for Z>20, due to 
increasing Coulomb barriers and decreasing 
Q-values. For heavier nuclei beyond Ca, the 
reaction pattern seems complicated, ffow¬ 
ever, temperatures approaching 10 9 K cause 
strong photodisintegrations for proton cap¬ 
ture Q-values smaller than 25kT&3 MeV, 
which is comparable to Q-values found for 
even-Z T z =-1 nuclei. [25-30 kT>Q is a rough 
criterion to find out whether photodisinte¬ 
grations can counterbalance capture reactions 
and lead to an equilibrium situation.] Thus, 
reactions breaking out of the cycles discussed 
above (up to the proton drip line) come into 
a (p, 7 )-( 7 , p)-equilibrium along isotonic lines 
with equal N. Under such circumstances the 
nuclear connection boils down to the nec¬ 
essary knowledge of nuclear masses and f3- 
decay half-lives, similar to the r-process ex¬ 
periencing an (n, 7 )-( 7 , n)-equilibrium in iso¬ 
topic chains (see |Cowan, Thielemann, fe Tru- 


rar f 199 1|; [Kratz ct al. 1993| ; |Meyer 1994|) . 

The general situation is explained in Fig¬ 
ure [| which contains all possible reactions 
from the original cycles over proton captures 
br eaking out of the cycles, the accompany¬ 


in g} beta, decays a.nrl (rvpi-reactions bridging 


waiting points with long decay half-lives. By 
comparing with Figure |2|, which also lists the 
Q-values involved, Figures [4a] and [44] give an 
indication for the minimum temperatures re¬ 
quired to permit statistical model calculations 
for proton and alpha-induced reactions. The 
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Fig. 3.— Recurrent structure in the rp- 

process network. 

reason is that a high density of compound nu¬ 
clear states is necessary for statistical model 
calculations to become reliable. The energy 
of the Gamow peak, where the integrand in 
a thermonuclear cross section integral (cross 
section times Boltzmann distribution) has a 
maximum, increases with increasing temper¬ 
ature. Because higher level densities are en¬ 
countered at higher excitation energies corre¬ 
sponding to higher temperatures of the Boltz¬ 
mann distribution, a critical temperature ex¬ 
ists beyond which such calculations are reli¬ 
able (for more details see [Rauschcr, Thiclm 
mann, fe Kratz] 1996). We can analyse the 
expected behavior: Reactions j, k, and l 
of Figure [3] have reaction Q-values less than 
~3 MeV and for conditions close to 10 9 K 
we have 25 kT^Q which permits a (p, 7 )- 
( 7 , p)-equilibrium. Thus, although the cross 
sections are highly uncertain and statistical 
model approaches are not valid (see Figure 
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Fig. 4a.— Applicability of the Statistical 

Model for (p, 7 )- and ( p , a)-reactions. 


Fig. 4b.— Applicability of the Statistical 
Model for ( 0 , 7 )- and (a,p)-reactions. 


ffa|), only the Q-values, i.e. masses, are needed 
for these nuclei. Reactions b, c, d , and e [(p, 7 ) 
and ( p , a) reactions] have Q-values in excess 
of 5 MeV and in most cases high enough 


lev |el densities to ensure the usage of statist 

tical model cross sections for the appropriate 
temperature regime. Reaction /, an (cp^- 
reaction, provides also a sufficiently high level 
density to apply statistical model cross sec¬ 
tions (see Figure [TE], for 18 Ne(a,p) see |Gorres 


& |Wiescher 1995|) . Reaction g is crucial for 
the break-out, a statistical model approach 
is clearly not permitted, and a (p, 7 )~( 7 , p)- 
equilibrium is not valid at the lower appropri¬ 
ate temperatures of about 3-4xl0 8 K. Thus, 
all reactions of type g ask strongly for experi¬ 
mental determination, possibly with radioac¬ 
tive ion beams. Reactions h, i and all other 
connecting /3-decays are required, preferably 
from experiments, otherwise from theoretical 
QRPA predictions. 

Beyond Ca and Ti, this scheme simpli¬ 
fies into one dominated by proton captures 
and beta-decays, with the proton captures in 


a (p, 7 )-( 7 , p)-equilibrium. The proton drip 
line develops a zig-zag shape and can separate 
two even -Z proton-stable nuclei on the same 
isotone by an unstable odd -Z nucleus. |Gorre: 


ct al.| (1995) have shown that such gaps can 


be bridged by two-proton capture reactions, 
proceeding through the unstable nucleus, sim¬ 
ilar to the 3a-reaction passing through 8 Be. 
This can extend the (p, 7 )^( 7 , p)-equilibrium 
beyond the drip line to proton-stable ’’penin¬ 
sulas”, thus avoiding long beta-decay half- 
lives. While alpha-induced reactions are not 
important for nuclei beyond Ca and Ti, small 
Q-values for A > 56 can permit ( 7 , a)-photo- 
disintegrations and cause a back flow, de¬ 
creasing the build-up of heavy nuclei. 

After having analized the reaction patterns 
in the previous discussion, we have now to 
implement these hirelings into an approxima¬ 
tion which can handle and describe all of the 
patterns. Considering the hydrogen and he¬ 
lium abundances as slowly varying, and as¬ 
suming a kind of local steady how by requir¬ 
ing that all reactions entering and leaving a 
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nucleus add up to a zero flux, we will devise 
a scheme which can adapt itself to the whole 
range of conditions of interest in hot hydro¬ 
gen burning. It can automatically cover situa¬ 
tions at low temperatures, characterized by a 
steady flow through all connected hot CNO- 
type cycles and branchings feeding the next 
cycle, to high temperature regimes where a 
(p, 7 )-( 7 , p)-equilibrium is established, while 
/3-decays or (a,p)-reactions feed the popula¬ 
tion of the next isotonic line. The details are 
discussed in the following section. 

3. APPROXIMATION SCHEME 



3.1. Reaction Rates and Quasi-Decay 
Constants 

The thermonuclear reaction rates used in 
our calculations have been discussed in detail 
by van Worrner et al. (|1994Q (c/. the ap¬ 
pendix therein) .Q The temporal change of the 
isotopic abundances Yj—Y. can be calculated 
by all depleting and producing reactions 




i 3,k 

£ Vi Y k Y,. 

j,k,l 


( 1 ) 


The first term in (I|) includes all one-particle 
reactions (decays or photodisintegrations) of 
nuclei j producing (i^j) or destroying ( i=j ) 
nucleus i. The aj are either decay con¬ 
stants or effective temperature dependent de¬ 
cay constants in the case of photodisintegra¬ 
tions. The second sum represents all two- 
particle reactions between nuclei j, k with 


lr The parameterized reaction rates can be found 
for instance in the database REACLIB located at 
|ittp://csa5. lbl.gov/ 'fchu/astro.html/. 


Fig. 5a.— rp-process network being used in 
the approximation. The structure in the box 
is eight times repeated up to 50 Mn. Bold ar¬ 
rows correspond to large, normal arrows to 
small decay constants. 

l oij ) k= l Cj,k pNA l (o'v)j t k- Similarly, the third 
term represents three particle processes with 
= l c jt k,ip 2 Nl l (av)j^i ([Fowler, Caughlan 


fc Zimmcrmannl 1967). The coefficients c in¬ 


clude the statistical factors for avoiding dou¬ 
ble counting of reactions of identical particles. 

van Worrner et al. ( |1994|) investigated 
numerically which reactions must be consid¬ 
ered in the rp-process. Apart from the tem¬ 
perature insensitive /3 + -decays, the only one- 
particle reactions which must be especially 
considered at high temperatures are photo¬ 
disintegrations. The two-body reactions dom¬ 
inating the rp-process are (p, 7 )-, (p, a)- and, 
(a, p)-reactions. Three body reactions can be 
neglected in explosive hydrogen burning with 
the exception of the 3a-reaction which plays 
a central role in this context. Taking into ac¬ 
count individual reactions, equation (]T|) can 
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Thus the system of nonlinear differential equa¬ 
tions can be written in a linearized form 

i’.=E*d y i. (4) 

3 

where x represents one of the following re¬ 
actions: ( 7 ,p), (p, 7 ), (a,p), ( p,a ), (a, 7 ), 
( 7 , a), ( 1 Q + v ). When we choose Y 0 =(Yi, ..., Y*, 
..., Y n ) to represent the initial abundance dis¬ 
tribution, the system of equations (||) can be 
written as 


Fig. 5b.— The proton capture sequences 
are connected to each other by (at,p)~ and 
(/? + , zz)-breakout reactions respectively. 

be rewritten as 

y < = eYY + EVY 

3 3 

+ E pNA^av^YjYp 

3 

+ T,pN A i (av) j , a Y j Y a . (2) 

3 

If the proton and helium abundances do not 
change greatly during relevant rp-process time 
scales, we are able to describe proton cap¬ 
tures and (a,p)-reactions with quasi-decay 
constants, provided that temperature and den¬ 
sity do not change greatly during the process 
time scale of interest: 



■= pN A i 

[° v )j, P Y P 

= const 

(3a) 

i\P’ a 

A 3 

■ = pN A i 

' av h, P Y p 

= const 

(3b) 

i \«>7 
X 3 

■ = pN A l 

W v )j, a Y a 

= const 

(3c) 

^7 

■■= pN A i 

'\ aV )j,a Y <x 

= const. 

(3d) 


Y = AY , A(N x Y) -matrix. (5) 

Such systems of linear differential equa¬ 
tions with constant coefficients can be solved 
by Jordan-normalizing matrix A. Its trans¬ 
formation 

J = SAS~ 1 ,Se GL n (R) ( 6 ) 

yields the TV-dimensional set of eigenfunctions 

u,j(t) ^ e^ p ,.( t ). (7) 

Here, A ” 1 represents eigenvalue i with degen¬ 
eracy m , Ui-j eigenvector j to eigenvalue i 
(1 <j <m) and Pij(t ) a polynomial of de¬ 
gree (j- 1). However, the transformations 
u' l . ] —Su l .j are in general difficult to survey 
and to handle. Furthermore, the elements of 
A become slowly varying functions when the 
abundances of hydrogen and helium change 
slightly over a certain process time scale. This 
is for instance the case, when we consider hot 
hydrogen burning. Only numeric simulations 
will be feasible for these two reasons. 

We will show in the following subsection 
how certain sequences of a thermonuclear re¬ 
action network can be regarded as quasi- 








abundance 




Fig. 6a.— Comparison of the approximated 
abundances of 53 Ni and 55 Ni with results from 
full network calculations. 

nuclei. In such sequences a steady flow ma¬ 
terializes and the ratios between the different 
abundances are constant and are only deter¬ 
mined by the values of the quasi-decay con¬ 
stants. A procedure of this kind provides us 
with a system of equations with a smaller di¬ 
mension, which can be solved numerically. 

3.2. The Approximation 

In section |], we discussed the burning cy¬ 
cles with the same structure as the hot CNO- 
cycle (see also Figures H>. It is well un¬ 
derstood how a steady state equilibrium is 
established in hot CNO-like reaction chains 
at temperatures T 9 ~ 0.2 and on time scales 
governed by the longest /3 + -decay in the cy¬ 
cle. Then the abundance Y hJ of nucleus ( i ) in 
a cycle (j) can be calculated from the equi¬ 
librium condition 

Y,j = 0 . ( 8 ) 

This condition can be used for simplifying 
the system of equations which determines the 
thermonuclear reaction flows. However, an 
approximation of this kind is no longer jus- 




Fig. 6b.— Same as Figure |6a] but for 64 Ge 
and 72 Kr. 

tided, when we consider time scales shorter 
than the half-lives of the /3 + -unstable nuclei or 
temperatures greater than T 9 =l where pho¬ 
todisintegrations become so important that 
a (p, 7)-(7, /^-equilibrium is attained rather 
than a steady flow (see Figure |3|). Then, 
proton-induced reactions can be associated 
with large quasi-decay constants and, in com¬ 
parison to /^-decays, small half-lives. Thus, 
all the abundances of isotonic nuclei which 
are connected by proton capture reactions 
can be considered in a quasistatic equilibrium 
on time scales exceeding the half-lives associ¬ 
ated with the quasi-decay constants. Such se¬ 
quences shall be now defined as quasi-nuclei. 

In the following paragraphs, the abundances 
of nuclei in proton capture sequences shall be 
calculated under the assumption of a quasi¬ 
static equilibrium. Consider the j th of m se¬ 
quences which consists of n —1 proton cap¬ 
tures between n nuclei (c/. Figure |H5|: Here, 
two proton capture sequences are shown with 
n=4 and n= 5). Each nucleus ( i ) in a sequence 
(j) is involved ink+2 reactions including pro¬ 
ton captures and photodisintegrations. Fur¬ 
thermore, we assume for simplicity that the 
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abundances of helium and hydrogen stay con¬ 
stant over the time scale of interest and that 
nucleus (i,j) is not fed by reactions from an- _ 
other chain j'. Using the equilibrium condi- '» 
tion (H), the temporal change of the abun- 

dance of a nucleus (i,j) can be written as “ 

0 ) 


V — \ p ’7 V — V” \ x V 

k 

+X]f 1J Y i+hj = 0 . (9) 


The terms summed over k are leaks out of 
chain (j) via nucleus (i,j) The extreme cases 
of flux equilibria, like steady flows on the one 
hand (A p ’ 7 A 7,p , A rj+ ) and (p, r y)-('y,p)- 

equilibria on the other hand (A p>7 , A 7,p 
A^ + ) can be discussed on the basis of this 
equation. Here, the general case shall be dis¬ 
cussed. 

For a chain of n nuclei, n —1 such equa¬ 
tions permit to relate the abundance U—ij to 
Yij starting with a “known” Y n j, the abun¬ 
dance of the last nucleus in a sequence. Thus 
we obtain for each sequence (j) an (n— 1)- 
dimensional system of linear equations which 
allows us to express each abundance Y tJ in 
terms of Y n j: 


YiJ = ■ ■ ■ K,j)Y„j . ( 10 ) 

The Vij depend only on the constant or slowly 
varying quasi-decay constants X t . r The abun¬ 
dance of a sequence can be treated as the 
abundance of a quasi-nucleus by Yf—Y^iY^j. 
When Yj is already determined at a point in 
time t=t s , each abundance Y t J can be written 
as 


Y id (t a ) = ^-Yj{t s ) , (11) 

v j 

with Vj—J2i Vij. Hence, when the abundances 
of the quasi-nuclei Yj are given at a time t, 



Fig. 7a.— Comparison of the energy genera¬ 
tion rate e for for nova-like conditions. 


the abundances Y tJ (t) in each sequence can 
be determined simultaneously. 

However, we also have to take care of in¬ 
coming fluxes into a sequence ( j ) not jet con¬ 
sidered in equation (|9|). Asa consequence, the 
coefficients v^j are no longer constant because 
they become explicitly abundance-dependent 
for the following reason: in contrast to the 
outgoing fluxes, the incoming fluxes into a se¬ 
quence j are not proportional to the abun¬ 
dances Yij of the nuclei in this sequence. 
Therefore, abundances can no longer be de¬ 
termined simultaneously. 

Incoming fluxes can nevertheless be consid¬ 
ered as long as the abundances do not change 
too fast during a time step. We can add 
to equation (|^) all incoming fluxes from nu¬ 
clei (l) of different chains (m) proportional to 
their abundances Yi m (t s _i), m^j (when no 
time t is indicated, the quantities are given at 
t=t s ), 

V — —V. \' V x I \7.p V. I \P>7 V. 

*i,j 2_^ A i,j,k ' A i+ 1-U+lj + A i-lj r i-l,j 

k 
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UJ 



Fig. 7b.— Comparison of the energy genera¬ 
tion e for nova-like conditions. 


+ E A£ m Yi,= 0 . ( 12 ) 

l,m 


This can be solved for leading to 


Yi -ij = 


TP^(^EA^-A?5rm,i) 
A-lj \ k J 


Vi — l j Y n 


\P,7 \,wXl,m{t s - 1) , (13) 

— l,jf l,m 
m^j 


or more generally 

^i,j Y n 'j Wi j(t s — i). (^4) 

When a slowly varying function u'(£ s _ 1 ) is 
given so that we can approximate Y n j by 


Y n,j(t s ) 


v'Ats- 1) 


, (15) 


we can 


extend w hJ in flTll) as follows 


Fig. 7c.— With this approximation, energy 
generation for nova-like conditions can be re¬ 
produced with an accuracy of 15-5 per cent. 


v- 

Yi,j = v i,jYn,j — TY u! i,jY n j , (16) 

1 3 

and finally obtain 

Yij = (17) 

/ v i t n 

v iJ = v i,3 - 77 Wij , (18) 

1 3 

with Vj—Y^ivY. Setting v' i j(t 0 =0)=Vij(t 0 =0) 
as initial values, the existence of the functions 
v'j is guaranteed. 

So far we have made two approximations, 
(a): including reactions within a chain and 
leaks out of a chain described by equation 
(|). In that case all abundances could be de¬ 
scribed simultaneously, (b): When also con¬ 
sidering incoming fluxes [see equation (|l^)], 
we become dependent on abundances in other 
chains. Approximating their abundances at 
t=t s by abundances at t'=t s - 1 permits to stay 
within the same approximation scheme. But 
we have to be aware that this is only valid 
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de/dt [erg g 



Fig. 8a.— Same as Figure [Fa| but for X-ray Fig. 8b.— Same as Figure [ft] but for X-ray 

burst conditions. burst conditions. 


for slowly varying abundances or equivalently 
small time steps. 

We obtained this approximation under the 
assumption that the abundances of hydrogen 
and helium stay constant. Y p and Y a are 
slowly varying functions in realistic environ¬ 
ments. We can adjust for the change of Y p 
and Y a and obtain their new values applying 
relation (|TTD 


Y„ = 


Y n = 


W/U,P I \Oi,p _ \P,a 
/ i.7 ' ^i,j 

X^(\P' a _ \ a ’P _L_ \^’ a 

-W)YiJ • ( 19b ) 


As a consequence, the quasi-decay con¬ 
stants must be recalculated after each time 
step, making use of the updated proton and 
alpha abundances considering them constant 
for one time step. 

Thus it is possible to express simultane¬ 
ously every abundance Y it] in a sequence with 


the sequence abundance Y y With an approx¬ 
imation scheme of this kind, we reduced our 
full network described by YTj li n (j) equations 
to a system of m equations. In the following 
we want to check the accuracy of such an ap¬ 
proximation. 

4. RESULTS 

The calculations were done with a reaction 
network which consists mainly of proton cap¬ 
ture sequences, each starting out from a T z =0 
nucleus. The number of nuclei in a sequence 
is determined by the condition that the quasi 
decay constant * +1 Af 7 of a nucleus (i) into 
nucleus (i+1) is larger than a certain A min . 
When this condition is fulfilled, nucleus (f+1) 
is considered to be in the sequence. For the 
following calculations A min was set to 10 3 s^ 1 
because we wanted to reproduce the energy 
generation of explosive hydrogen burning on 
time scales r=A” 1 > 10~ 2 —10 _1 s. As a con- 

min 

sequence, a sequence contains not more than 
5 nuclei for temperatures Tg<1.5 and densi¬ 
ties p<10 6 g cm -3 . 

In addition to the proton capture sequences, 
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Fig. 8 c.— Same as Figure [7c| but for X-ray 
burst conditions. 


T, = | nuclei had to be considered in the net¬ 
work. They are the most neutron-rich nu¬ 
clei of the hot CNO-type cycles, where after 
a proton capture, the alpha-threshold in the 
compound nucleus is lower than the proton- 
threshold, and (p, a)-reactions can occur which 
compete with the still possible (p, 7 )-transitions 
(see Figure |5t]). Since the (p, a)- and (p, 7 )- 
reactions occur on the same time scale, these 
T z =\ nuclei can not affiliated to a proton cap¬ 
ture chain but have to be considered sepa¬ 
rately. 

Between 12 C and 50 Mn, the network can 
be composed of proton capture reaction se¬ 
quences and T z =\ nuclei with one exception. 
The transition occurs already at compound 
nucleus 19 Ne in the Ne isotopic chain. A 
(p, a)-reaction is possible on the T z =0 nucleus 
18 F and a proton capture reaction sequence 
starts out from 19 Ne with isospin T z —— |. 

The low alpha-threshold of 19 Ne, which 
caused an exception of the rule and prevented 
an OFNe-cycle, permits on the other hand the 
15 0 (a, 7 )-reaction and triggers the break-out 
from the CNO-cycles. The network described 


P=P 0 (T/T/A T=T Q e- 
r = 120 s, p =10 4 g/cm 3 , T =0.25xl0 9 K 

' r q =>/ ’ 0 



Fig. 9.— Same as Figure [Fa] but for varying 
density and temperature starting with novae 
peak conditions. The approximation is valid 
as long as the temperature does not fall con¬ 
siderably below Xg = 0 . 1 . 


so far is shown in the Figure |5a| and [X]. 

At temperatures exceeding T 9 = 0.8 and p > 
10 6 g cm -3 , the network must be completed 
with the 3a-reaction. Furthermore, it is very 
likely that fluxes up to ' 2 Kr may influence 
the energy generation. I 11 the region between 
51 Mn and ' 2 Kr, the /3 + -reactions of the for 
the rp-process relevant nuclei are, with the 
exception of 55 Ni, 64 Ge and ' 2 Kr, consider¬ 
ably smaller than the typical time scales of 
novae and X-ray bursts. Consequently, only 
the above mentioned nuclei were considered 
explicitly in the approximation scheme. The 
approximated abundances of 53 Ni and these 
three waiting point nuclei are compared with 
full network calculation results in Figure | 6 a| 
and BE. 


The nuclei used in the approximation scheme 
are listed in Table 1 according to their af¬ 
filiation to a certain proton capture reaction 
chain. Nuclei which had to be considered sep¬ 
arately (p, 4 He, T z —\ nuclei, 55 Ni, 64 Ge and 
, 2 Kr) are given an own sequence number. 
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P=P 0 ( T /' T o) 4/3 . T=T 0 e-^ 

r = 1.2 s, p =10 6 g/ctn 3 , T =1.5xl0 9 K 
' 0 te/ ’ 0 



Fig. 10.— Same as Figure [| but for X-ray 
burst peak condition. 


temperatures and densities, i. e. 10 4 <p< 10 6 , 
0.2<T g <1.5 confirmed the above mentioned 
efficiency and accuracy of the approximation 
scheme. 

Such an approximation would not be very 
useful if it was only valid for constant tem¬ 
peratures and densities. In order to check its 
accuracy for varying temperatures and den¬ 
sities we assumed after peak conditions were 
obtained in the ignition an exponentially de¬ 
creasing temperature and an adiabatically ex¬ 
panding accretion layer with a radiation dom¬ 
inated equation of state (7 = |): 


The waiting point abundances obtained with 
this approximation are compared with full 
network calculation results in Figure ^a] and 

B 


T = T 0 e * 

( 20 a) 

/T\i 


p = P0 {tJ ■ 

( 20 b) 


Approximating the network according to 
section pTT2| . we calculated the energy gener¬ 
ation rate e and the integrated energy gen¬ 
eration e for nova-like conditions (Tg=0.2, 
p=10 5 g cm -3 ) and X-ray burst peak condi¬ 
tions (Tg—1.5, p=10 6 g cm -3 ) and compared 
the results to full size network calculations 
which take into account each nucleus sepa¬ 
rately. In each case we assumed a solar abun¬ 
dance distribution for the accreted matter. 
The results for small temperatures are shown 
in Figure [7a] [A] and those for X-ray burst con¬ 
ditions in Figure [ 8 a| [ 8 c|. Considering the time 
scales important for novae ( 100-1000 s) and 
X-ray bursts (1-10 s), we see that the approx¬ 
imation scheme devised in section 3 agrees 
always to a high degree for the astrophysi- 
cal applications considered here. The overall 
results show only deviations smaller than 15 
per cent or even less (down to 5 per cent) 
at a gain in computational speed by a fac¬ 
tor of 15. Several calculations at intermediate 


Using typical peak temperatures T 0 , densi¬ 
ties p 0 and time scales r for X-ray bursts and 
novae, we again calculated the energy genera¬ 
tion (Figure || and [10]). The results show that 
the approximation is valid for changing tem¬ 
peratures and densities as long as they are 
in the range of 10 4 <p<10 6 , 0.1<T 9 <1.5. No 
significant change in computational speed was 
observed. 


The differences at time scales smaller than 
f=lCT 2 s are due to the very fast reaction 
flows within the proton capture reaction chains: 
They occur in full scale network calculations 
until equilibrium is reached according to equa¬ 
tion (||). 


Table 1: 

Proton Capture Sequences 1 
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Sequence Number 

Nuclei 

Sequence Number 

Nuclei 

Sequence 1 .. 


P 

Sequence 6 . 

^Ne" 

Sequence 2 .. 


4 He 


20 Na 

Sequence 3 .. 


12 C 


21 Mg 



13 N 

Sequence 7 . 

16 0 



14 0 


17 f 

Sequence 4 .. 


14 N 


18 Ne 



15 0 

Sequence 8 . 

20 Ne 

Sequence 5 .. 


15 N 


21 Na 





22 Mg 

Table 1 - 

- Continued 


23 A1 


24 Si 

Sequence 9 . 22 Na 

23 Mg 
24 A1 
25 Si 


Sequence 10 . 23 Na 

Sequence 11 . 24 Mg 


25 A1 

26 Si 

27p 

28 S 


Sequence 12 . 26 A1 

27 Si 

28 p 

29 S 

Sequence 13 . 2 'A1 

Sequence 14 . 28 Si 

29 p 

30 q 

31 a 


32 Ar 


Sequence 15 . 30 P 

31 S 

32 Cl 

33 Ar 

Sequence 16 . 31 P 
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Table 1 — Continued 


Table 1 — Continued 


Sequence Number 

Nuclei 

Sequence Number 

Nuclei 

Sequence 17 . 

32 S 

Sequence 27 . 

45 V 


33 C1 


47 Cr 


34 Ar 


£ 

00 


35 K 


49 Fe 


36 Ca 

Sequence 28 . 

47 V 

Sequence 18 . 

34 C1 

Sequence 29 . 

48 Cr 


35 Ar 


49 Mn 


36 K 


50 Fe 


37 Ca 


51 Co 

Sequence 19 . 

35 d 


52 Ni 

Sequence 20 . 

36 Ar 

Sequence 30 . 

50 Mn 


37 K 


51 Fe 


38 Ca 


52 Co 


39 Sc 


53 Ni 


4o Ti 

Sequence 31 . 

51 Mn 

Sequence 21 . 

38 K 

Sequence 32 . 

55 Ni 


39 Ca 

Sequence 33 . 

64 Ge 


40 Sc 

Sequence 34 . 

72 Kr 


41 Ti 

^The nuclei used in the approximation scheme 

Sequence 22 . 

39 K 

are listed according to their 

affiliation to a 

Sequence 23 . 

40 Ca 

certain proton capture reaction chain. Nuclei 


41 Sc 

which had to be considered separately (p, 4 He, 


42 Ti 

T z = i nuclei, 55 Ni, 64 Ge and 

72 Kr) are given 


43 V 

an own sequence number. 



44 Cr 




Sequence 24 . 42 Sc 5. CONCLUSION 

43 Ti 

44y The main motivation for the present in- 

450. vestigation was to find a fast and efficient 

Sequence 25 43 Sc approximation scheme which permits to pre- 

Sequence 26 . 44 Ti diet the energy generation in explosive hydro- 

45y gen burning for a large range of conditions, 

46Q r i.e. temperatures (0.2<T 9 <1.5) and densities 


47y[ n (10 4 g cm -3 <p< 10 6 g cm -3 ) respectively, for 

48pg applications in hydro calculations which can¬ 

not afford full nuclear network. This goal has 
been achieved. The energy generation of full 
network calculations (150 nuclei) could be re¬ 
produced with a high accuracy, resulting in 
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p= 10 6 g cm 3 , T =°1.5 



Fig. 11.— Comparison of the 33 Ar- 

abundances calculated by a full network cal¬ 
culation and the approximation scheme. 

deviations of 5 to a maximum of 15 per cent 
while gaining a factor of about 15 in compu¬ 
tational speed. 

The approximation scheme discussed in 
this paper is therefore well suited for realistic 
hydro calculations of novae or X-ray bursts 
and other possible sites of explosive hydrogen 
burning. 

The additional advantage of this method is 
that it also permits to predict isotopic abun¬ 
dances with a similar accuracy and can thus 
even replace full network calculations for this 
purpose (see Figure 

We were also able to determine key nu¬ 
clear properties, which directly enter the pre¬ 
cision of calculations for explosive hydrogen 
burning, its energy generation and abundance 
determination. We have identified the even- 
Z T z =-l/2 nuclei like 23 Mg, 27 Si, 31 S, 35 Ar 
and 39 Ca, as essential targets for proton cap¬ 
tures which are in competition with / 1 -decays. 
As their small reaction Q-values (less than 
2 MeV) do not permit the application of sta¬ 


tistical model cross sections, they have to 
be determined experimentally. Thus, these 
results are an ideal guidance for future ex¬ 
perimental investigations with radioactive ion 
beams. Nuclei closer to stability permit the 
application of statistical model cross sections. 
Nuclei more proton-rich than T z —~ 1 are only 
populated for conditions when a (p, 7 )^( 7 , p)- 
equilibrium is established. 

Thus, only their masses, spins and half- 
lives enter the calculation. Connecting (ab¬ 
reactions have high enough Q-values, even 
close to the drip-line, that the level densi¬ 
ties are sufficient for statistical model cross 
sections (this is true at least up to Ca and 
Ti, where alpha-induced reactions can play 
an important role). Therefore, the remain¬ 
ing experimental properties which should be 
determined and enter the calculations in a 
crucial way are /3 + -decay half-lives and nu¬ 
clear masses. The latter control the abun¬ 
dances in a (p, 7 )-( 7 ,p)-equilibrium. Beyond 
Se, the details of the proton drip-line are ac¬ 
tually not that well known yet, but can influ¬ 
ence the endpoint of the rp-process. First cal¬ 
culations which include 2 p-capture reactions 
(Schatz et al. 1996) indicate that it seems 
possible to produce nuclei with A=90-100 in 
the rp-process under X-ray burst conditions, 
depending on the choice of mass formulae and 
half-life predictions (similar to r-process stud¬ 
ies). Especially a region of strong deformation 
and small a-capture Q-values around A=80 
seems to be predicted with large variations 
among different mass models, and the even- 
even N=Z nuclei between A =68 and 100 play 
a dominant role. 
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